Comparison between the sexes (Male vs Female) for the ages and microbiota

Table of content

  1. Differentially expressed genes
  2. Stats about the DEG
  3. DEG with significant p-value and fold change
    1. Log2FC
    2. Z-score
  4. Co-expression (WGCNA))
    1. Z-score in modules
    2. Genes in modules
  5. Enrichment analysis
    1. GO analysis
    2. KEGG pathways

Generated from a Jupyter Notebook - Sources

Loads

Libraries and functions

Warning message in is.na(x[[i]]):
“is.na() applied to non-(list or vector) of type 'environment'”Warning message in rsqlite_fetch(res@ptr, n = n):
“Don't need to call dbFetch() for statements, only for queries”
==========================================================================
*
*  Package WGCNA 1.63 loaded.
*
*    Important note: It appears that your system supports multi-threading,
*    but it is not enabled within WGCNA in R. 
*    To allow multi-threading within WGCNA with all available cores, use 
*
*          allowWGCNAThreads()
*
*    within R. Use disableWGCNAThreads() to disable threading if necessary.
*    Alternatively, set the following environment variable on your system:
*
*          ALLOW_WGCNA_THREADS=<number_of_processors>
*
*    for example 
*
*          ALLOW_WGCNA_THREADS=4
*
*    To set the environment variable in linux bash shell, type 
*
*           export ALLOW_WGCNA_THREADS=4
*
*     before running R. Other operating systems or shells will
*     have a similar command to achieve the same aim.
*
==========================================================================


Allowing multi-threading with up to 4 threads.
[1] "preparing gene to GO mapping data..."
[1] "preparing IC data..."
[1] "preparing gene to GO mapping data..."
[1] "preparing IC data..."
[1] "preparing gene to GO mapping data..."
[1] "preparing IC data..."

Data

Differentially expressed genes

Extract DEG between Male and Female for the different ages and microbiota combinations

  • Threshold for adjusted p-value: 0.05
  • Threshold for adjusted significant fold change: 1.5

Table with the factors / constrasts

InfoInterceptMale vs FemaleGF vs SPFMiddle-aged vs YoungOld vs YoungMale & Middle-agedMale & OldMale & GFGF & Middle-agedGF & Old
Male vs Female (SPF, Young) 0 1 0 0 0 0 0 0 0 0
Male vs Female (GF, Young) 0 1 0 0 0 0 0 1 0 0
Male vs Female (SPF, Middle-aged)0 1 0 0 0 1 0 0 0 0
Male vs Female (GF, Middle-aged) 0 1 0 0 0 1 0 1 0 0
Male vs Female (SPF, Old) 0 1 0 0 0 0 1 0 0 0
Male vs Female (GF, Old) 0 1 0 0 0 0 1 1 0 0

Extract the log2FC of the DEG

Stats

Using type as id variables
Male vs Female (SPF, Young)Male vs Female (GF, Young)Male vs Female (SPF, Middle-aged)Male vs Female (GF, Middle-aged)Male vs Female (SPF, Old)Male vs Female (GF, Old)
All DEG (Wald padj < 0.05)2414303735673593261 227
All over-expressed genes (Wald padj < 0.05 & FC > 0)1338160018221802141 112
All under-expressed genes (Wald padj < 0.05 & FC < 0)1076143717451791120 115
DEG (Wald padj < 0.05 & abs(FC) >= 1.5)1291142916531642116 143
Over-expressed genes (Wald padj < 0.05 & FC >= 1.5)1085118412711233 65 64
Under-expressed genes (Wald padj < 0.05 & FC <= -1.5) 206 245 382 409 51 79

All DEG (Wald padj < 0.05)

DEG (Wald padj < 0.05 & abs(FC) > 1.5)

pdf: 2

DEG (Wald padj < 0.05 & abs(FC) > 1.5)

Log2FC

Z-score

Column order: age - microbiota - sex

Column order: microbiota - age - sex

Co-expression (WGCNA)

Z-score in modules

Column order: age - microbiota - sex

Column order: microbiota - age - sex

Genes in modules

Enrichment analysis

Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”Warning message in pcls(G):
“initial point very close to some inequality constraints”

GO analysis

Warning message in stack.default(getgo(l$sign_fc_deg$genes, "mm10", "geneSymbol")):
“non-vector elements will be ignored”

Biological process

Dot-plot with the 10 most significant p-values for the different comparison

Using category as id variables
Using category, type as id variables
Warning message:
“Column `variable` joining factors with different levels, coercing to character vector”

Cellular components

Dot-plot with the 10 most significant p-values for the different comparison

Using category as id variables
Using category, type as id variables
Warning message:
“Column `variable` joining factors with different levels, coercing to character vector”

Molecular functions

Dot-plot with the 10 most significant p-values for the different comparison

Using category as id variables
Using category, type as id variables
Warning message:
“Column `variable` joining factors with different levels, coercing to character vector”

GO networks

The edge colors in the tree represent the relationship between two nodes. - green: positively regulates

  • red: negatively regulates
  • black: regulates
  • blue: is a
  • light blue: part of

GO Trees at "../results/dge/type-effect/type_gender_age/go/"

[1] "Male vs Female (SPF, Young)"
Warning message:
“Column `values` has different attributes on LHS and RHS of join”Warning message:
“Column `values` has different attributes on LHS and RHS of join”
[1] "Male vs Female (SPF, Middle-aged)"
Warning message:
“Column `values` has different attributes on LHS and RHS of join”Warning message:
“Column `values` has different attributes on LHS and RHS of join”
[1] "Male vs Female (SPF, Old)"
Warning message:
“Column `values` has different attributes on LHS and RHS of join”Warning message:
“Column `values` has different attributes on LHS and RHS of join”
[1] "Male vs Female (GF, Young)"
Warning message:
“Column `values` has different attributes on LHS and RHS of join”Warning message:
“Column `values` has different attributes on LHS and RHS of join”
[1] "Male vs Female (GF, Middle-aged)"
Warning message:
“Column `values` has different attributes on LHS and RHS of join”Warning message:
“Column `values` has different attributes on LHS and RHS of join”
[1] "Male vs Female (GF, Old)"
Warning message:
“Column `values` has different attributes on LHS and RHS of join”Warning message:
“Column `values` has different attributes on LHS and RHS of join”

KEGG pathways

Warning message in stack.default(getgo(l$sign_fc_deg$genes, "mm10", "geneSymbol", :
“non-vector elements will be ignored”
Error in `$<-.data.frame`(`*tmp*`, labels, value = c("", "", "", "", "", : replacement has 34 rows, data has 48
Traceback:

1. plot_kegg_pathways(deg$KEGG$over$category, deg$sign_fc_deg, paste("../results/dge/", 
 .     dir_path, "/kegg/over_repr_kegg/", sep = ""))
2. suppressMessages(pathview(gene.data = fc_deg, pathway.id = cat, 
 .     species = "Mus musculus", gene.idtype = "Symbol"))
3. withCallingHandlers(expr, message = function(c) invokeRestart("muffleMessage"))
4. pathview(gene.data = fc_deg, pathway.id = cat, species = "Mus musculus", 
 .     gene.idtype = "Symbol")
5. `$<-`(`*tmp*`, labels, value = c("", "", "", "", "", "", "", 
 . "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", 
 . "", "", "", "", "", "", "", "", "", "", ""))
6. `$<-.data.frame`(`*tmp*`, labels, value = c("", "", "", "", "", 
 . "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", 
 . "", "", "", "", "", "", "", "", "", "", "", "", ""))
7. stop(sprintf(ngettext(N, "replacement has %d row, data has %d", 
 .     "replacement has %d rows, data has %d"), N, nrows), domain = NA)

Pathway graphs available at ../results/dge/gender-effect/gender_type_age/over_repr_kegg/

Pathway graphs available at ../results/dge/gender-effect/gender_type_age/under_repr_kegg/

Citations